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Recent cosmological data have provided evidence for a "dark" relativistic background at high 
statistical significance. Parameterized in terms of the number of relativistic degrees of freedom A'^off, 
however, the current data seems to indicate a higher value than the one expected in the standard 
scenario based on three active neutrinos. This dark radiation component can be characterized 
not only by its abundance but also by its clustering properties, as its effective sound speed and 
its viscosity parameter. It is therefore crucial to study the correlations among the dark radiation 
properties and key cosmological parameters, as the dark energy equation of state or the running 
of the scalar spectral index, with current and future CMB data. We find that dark radiation with 
viscosity parameters different from their standard values may be misinterpreted as an evolving dark 
energy component or as a running spectral index in the power spectrum of primordial fiuctuations. 

PACS numbers: 98.80.-k 95.85.Sz, 98.70.Vc, 98.80.Cq 
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I. INTRODUCTION 

The cosmological abundance of relativistic and " dark" 
particles as active or sterile neutrinos or any other light 
particle (as axions, for instance) is usually parameterized 
in terms of the number of relativistic degrees of freedom 
iVeff, being N^s = 3.046 the standard value for three 
(active) neutrino species. 

Several recent cosmological data analyses have con- 
strained the value of TVeff with increasing accuracy (see, 
for instance, Refs. [IHl]) providing a clear evidence for 
a dark relativistic background at high statistical signif- 
icance. The latest measurements of Cosmic Microwave 
Background (CMB) anisotropics at arc-minute angular 
scales from the South Pole Telescope (SPT) Q and the 
Atacama Cosmology Telescope (ACT) [l^, when com- 
bined with other cosmological data sets, yield the con- 
straint TVoff = 4.08^2,68 at 95% confidence level (CL) [p 
(see also [l^ and [l3|), showing evidence for dark radia- 
tion (i.e. NcS > 0) at more than seven standard devia- 
tions, suggesting values higher than those expected in the 
standard scenario. This hint, that clearly must be tested 
by future CMB measurements as those expected, for ex- 
ample, by the Planck satellite [l^, is interesting since an 
extra relativistic component in the standard three active 
neutrino model could be explained by a sterile neutrino. 
Models with one additional ~ 1 eV massive sterile neu- 
trino, i.e. the so-called (3-1-1) models, were introduced 
to explain LSND short baseline antineutrino data [l^ by 
means of neutrino oscillations. Up to date cosmologi- 
cal constraints on the number of massive sterile neutrino 
species have been presented in Refs. [1, [l6l - [l8| . How- 
ever, a larger value for A'eff could also arise from different 
physics, related to axions [l^, decay of non- relativistic 



matter [20|, gravity waves [21|, extra dimensions 
dark energy |23j or asymmetric dark matter models 
The constraints on iVoff are also affected by possible in- 
teractions of the dark radiation c omp onent with the dark 
matter component, see Refs. |24| - [2a |. 

Information on the dark relativistic background, how- 
ever, can be obtained not only from its effects on the 
expansion rate of the universe but also from its cluster- 
ing properties. For example, if dark radiation is made 
of massless neutrinos it should behave as relativistic par- 
ticles also from the point of view of perturbation the- 
ory. Following the definition introduced by Hu [13], this 
means that the dark radiation component should have an 
effective sound speed c^j^ and a viscosity parameter c^jg 
such that Cgjj = c^jg = 1/3. These perturbation param- 
eters can be constrained through measurements of the 
CMB anisotropics since dark radiation is coupled trough 
gravity with all the remaining components psj . This 
clearly opens a new window for testing the dark radia- 
tion component, since, for example, a smaller value for 



^ There is no fundamental symmetry in nature forcing a definite 
number of right-handed (sterile) neutrino species, as those are 
allowed in the Standard Model fermion content. 



c^jg could indicate possible non standard interactions (see 
e.g. [111). A value of c^jg different from zero, as expected 
in the standard scenario, has been detected in [30| and 
confirmed in subsequent papers [ll|, [3l| . More recently, 
a value of the effective sound speed c^^ smaller than the 
standard value of 1/3 has been claimed in Ref. [Tsj . 

From this discussion is clear that current cosmologi- 
cal data analyses are sensitive to dark radiation proper- 
ties and that the latest constraints on these parameters 
are showing interesting deviations from their expected 
standard values. It is therefore timely to investigate the 
impact of non-standard dark radiation properties in the 
determination of cosmological parameters related to dif- 
ferent sectors as dark energy or infiation. 

In this paper we derive bounds from the most recent 
cosmological data on the dark radiation parameters A^eff, 
c^ff, and c^jg relaxing the usual assumption of a ACDM 
cosmology, analysing the correlations among the cluster- 
ing parameters c^g and c^jg and other key cosmological 
parameters, as the dark energy equation of state or the 
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scalar spectral index. We also study these correlations 
in light of future Planck and COrE CMB data. We shall 
generate mock CMB data for a cosmology with dark radi- 
ation perturbation parameters different from their stan- 
dard values and then we fit these simulated data via 
the usual MCMC analysis to an extended non minimal 
cosmology with standard dark radiation parameters, but 
varying both the constant and the time varying dark en- 
ergy equation of state, or the scalar spectral index and 
its running. The paper is organized as follows. Section 
im describes the details of the analysis carried out here, 
including the cosmological parameters and data sets used 
in the analyses. In section Hill the different cosmological 
scenarios considered here are analyzed and the most im- 
portant degeneracies among the parameters are carefully 
explored. Section [IVI summarizes our main conclusions. 



II. ANALYSIS METHOD 

We perform our analyses considering three different 
scenarios: we first analyze the Wilkinson Microwave 
Anisotropy Probe (WMAP) data together with the lu- 
minous red galaxy clustering results from the Sloan Dig- 
ital Sky Survey II (SDSSII) ^ and with a prior on the 
Hubble constant Hq from the Hubble Space Telescope 
(HST) [Hi, referring to it as the "Baseline model". Then 
we add to these data sets SPT data [§| and we will refer to 
this case as the "BaselineSPT model" . We will explore as 
well the impact of Supernovae la (SNIa) luminosity dis- 
tance measurements from the Union 2 compilation [3J| in 
constraining the dark radiation parameters, and we will 
refer to this case in the following as the "BaselineSPT- 
SNIa model" . 

We have modified the Boltzmann CAMB code [s^ in- 
corporating the two extra dark radiation perturbation 
parameters c^g and c^jg and we have extracted the cos- 
mological parameters from current data using a Monte 
Carlo Markov Chain (MCMC) analysis based on the pub- 
licly available MCMC package CosmoMC [s^. We sample 
the following six-dimensional standard parameters: the 
baryon and cold dark matter densities (wb = fli,h^ and 
u!c = flch'^), the ratio between the sound horizon and the 
angular diameter distance at decoupling O^, the optical 
depth T, the scalar spectral index Us, and the ampli- 
tude of the primordial spectrum . We consider purely 
adiabatic initial conditions and we impose spatial flat- 
ness. We also include the effective number of relativistic 
de grees of freedom N^g, the effective sound speed c^^^ 
and the viscosity parameter c^j^. Table U shows the flat 
priors considered on the different cosmological parame- 
ters. Finally, we generate a mock data set for the ongo- 
ing Planck CMB experiment, with c^jj, different from its 
standard value and with w = —1 and Hs ~ 0.96. Then 
we fit these mock data using a MCMC analysis to dif- 
ferent extensions of the minimal cosmological model in 
which the dark radiation is standard. The three possi- 
ble extensions we consider are: (a) a. ACDM model with 



a running spectral index n^un, (b) the luCDM model in 
which we include the possibility of a dark energy equa- 
tion of state parameter w different from —1, and (c) the 
u'(a)CDM model in which we assume an equation of state 
evolving with redshift. The reconstructed values of the 
dark energy equation of state and of the running spectral 
index will be, in general, different from the values used 
in the mocks and, in the case of the dark energy equa- 
tion of state w, different from the value expected within 
the ACDM model. We shall also explore the impact of 
future CMB data from the COrE mission [s^] , perform- 
ing an equivalent forecast to the one we present here for 
Planck. 



Parameter 


Prior 




0.005 - 


0.1 




0.01 - 


> 0.99 




0.5 - 


> 10 


T 


0.01 - 


^ 0.8 


Us 


0.5 - 


> 1.5 




2.7 - 


4 


iVeff 


- 


> 9 


Wis 


- 


> 1 


2 
Coff 


- 


> 1 




-0.2 - 


0.1 




-2 - 





■Wa 


-1 - 


^ 1 



TABLE I. Flat priors for the cosmological parameters consid- 
ered here. 



III. RESULTS 

A. Baseline models 

We consider a cosmological model described by the fol- 
lowing set of parameters: 

{ujb, wc, Qs, T, w, Us, log[10^°yls], Ncs, c^is, c^g}, 

Notice from the results in the first column of Tab. [H] 
that in the "Baseline model", (i.e. the one with 
WMAP7-hSDSSII-hi?o data), the preferred value for the 
effective number of relativistic degrees of freedom is 
A^off = 5-82lg'g4, considerably higher than the standard 
model prediction A^eff = 3.04. The addition of SPT data 
in the "BaselineSPT model" decreases the preferred value 
of Aoff, making it closer to its canonical value. Finally, 
the addition of SNIa data in the "BaselineSPT-SNIa 
model" brings the value of Ncs even closer to 3.04, how- 
ever, it is still ~ 2tT away from the former value. These 
results seem to agree with the excess of radiation claimed 
in the literature in previous analyses [11-0, d, [U . Re- 
garding the dark radiation perturbation parameters c^jg 
and c^ff , their preferred values are close to their standard 
values of 1/3, being c^g much better constrained from 
current data than A^. 
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There exists a degeneracy between c^g and Ug and 
between c^j^ and rig, see the top and middle panels of 
Fig. [TJ These degeneracies are related to the fact that 
if Us increases (decreases), the power at low multipoles 
decreases (increases) , while the power at high £ increases 
(decreases) to a lesser extent. This effect could be com- 
pensated by an increase (decrease) in the viscosity pa- 
rameter c^jg, that leads to a decrease (increase) of the 
power at all scales. Concerning the degeneracy between 
c^g and rig, it is mainly related to the degeneracy be- 
tween Cgjj and c^jg (see the lower panel of Fig. [T]). There 
also exist degeneracies between w and c^jg and between 
w and iVeff, see Fig. [2l As we shall explain in the fol- 
lowing section, a value of w > — 1 shifts the positions of 
the CMB acoustic peaks to lower multipoles £; this effect 
could be compensated by a decrease of c^jg or by an in- 
crease of Ncs- The degeneracy between w and c^jg gets 
alleviated when information on high i multipoles from 
SPT is considered in the analysis. 

A change on the scalar spectral index Ug can also be 
compensated by a change in N^s, see the upper panel of 
Fig. El This degeneracy only affects the amplitude of the 
CMB peaks: a higher A'off will reduce the amplitude of 
the CMB peaks at £ > 200 due to a higher Silk damping, 
which in turn is due to the increased expansion rate |12| . 



B. Future Planck and COrE CMB data analysis 
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In the following we shall present the reconstructed val- 
ues of Us, rirun, wq and Wa which will result from a 
fit of Planck and COrE mock data (generated with non 
standard values for the dark radiation perturbation pa- 
rameters, c^jg = 0.1) to a cosmology with a standard 
value for the dark radiation parameter c^jg = 1/3 but 
with a running spectral index or a time varying dark en- 
ergy component. We do not consider here c^g 7^ 1/3 due 
to the tighter current bounds on this parameter, when 
compared to the current constraints on c^. . 



ACDM + riyun 




0.6 



For this scenario we consider the following set of pa- 
rameters: 

{uJt, UJc, ©s, T, Ug, \og[10^°As], nrun} • 

In general, the spectrum of the scalar perturbations is not 
exactly a power law but it varies with scale. Therefore 
one must consider the scale dependent running of the 
spectral index rirun = dus/dlnk. Following (38|, the 
power spectrum for the scalar perturbations reads 

/ 7 \ Us A- ^^{k/ ko){dn/dln k) ~\- ••• 

P{k) ^ Agfc"W oc (jA 
being = 0.05 Mpc~^ the pivot scale. The correlation 



FIG. 1. The top, middle and bottom panels depict the 68% 



'-Vis 



Us and 



„2 
Cefl 



and 95% CL bounds in the c^ff 
planes, respectively. The blue, red and green regions refer 
to the Baseline, BaselineSPT and BaselineSPT-SNIa models, 
respectively. 



between and nrun is shown in Fig.U] As stated in j39| . 
the parameter that is constrained by cosmological data is 
the effective spectral index n' = ns + \\i{k/kQ){dn/ dhik). 
This is the reason for the circular allowed regions in the 
IT'S — ?^iun plane. The first and second columns of Tab. lIIll 
show that, if a cosmology with rij-un — but with non 
standard dark radiation perturbation parameters (c^ 
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Baseline 
model 


BaselineSPT 
model 


BaselineSPT-SNIa 
model 


w 


—0.76 ± 0.15 


—0.85 ± 0.12 


—0.85 ± 0.12 


AT 


r on+O.60+2.71 
■^■"^-0.84-2.12 


- qo+0.27+1.07 
^••^"-0.31-0.98 


. r,Q+0.26+1.05 
^■^='-0.31-0.96 


c2. 

VIS 


rji +0.10 + 0.21 

'^•^J^-o.io-o.is 


n 0/1+0. 032+0. 17 
^•^^-0.052-0.13 


n oc:+0.03+0.42 
^•^•^-0.06-0.13 


en 


„ OI.+0.01+0.05 
'^•'-''^-0.02-0.05 


n .10+0.006+0.024 
'^••^'^-0.007-0.024 


r, O.1+0.01+0.02 

'^•■'•^-0. 01-0. 03 


Us 


0.976 ± 0.026 


0.982 ±0.024 


0.980 ± 0.024 



TABLE II. Constraints on the cosmological parameters for the three Baseline models described in the text. We report the 68% 
and 95% CL limits for the dark radiation parameters, and the mean and the standard deviation of the posterior distribution 
for the other cosmological parameters. 



ACDM + Urun (Planck) ACDM + Umn (COrE) wCDM (Planck) wCDM (COrE) «)(a)CDM (Planck) w(a)CDM (COrE) 



w -1 -1 -0.70 ±0.05 -0.63 ±0.05 — — 

N^s 3.04 3.04 3.04 3.04 3.04 3.04 

Us 1.002 ±0.004 0.999 ±0.002 1.007 ± 0.004 1.004 ± 0.002 1.007 ±0.004 1.007 ± 0.002 

n^un -0.035 ±0.005 -0.038 ± 0.003 

wo — — — — -1.19 ±0.10 -0.99 ±0.05 

Wa — — — — 0.77 ±0.23 0.88 ±0.06 



TABLE III. Constraints on the cosmological parameters for each of the Plank and COrE mock data sets described in the text. 
We report the mean and the standard deviation of the posterior distribution. We have set c^g — 1/3 and cjjg = 0.1 in the 
mock data sets used as fiducial models. Then, we have fitted these data to a model with canonical values for the dark radiation 
perturbation parameters, i.e. Cgft = 1/3 and Cvis = 1/3. 



0.1) is fitted to a cosmology with standard dark radiation 
parameters but with nrun 7^ 0, the reconstructed value of 
the running spectral index will differ from zero at a high 
statistical significance. 

Finally, for the case of the simulated cosmology here 
with c^jj, < 1/3, the reconstructed value of ris is consis- 
tent with a Harrison- Zel'dovich scale invariant primordial 
power spectrum within one sigma. Setting the proper- 
ties of dark radiation is therefore mandatory since it is 
highly correlated with the spectral index of the spectrum 
of primordial perturbations, key to distinguish among the 
different inflationary models. 



C. toCDM 

Here we consider a cosmological model including a dark 
energy fluid characterized by a constant equation of state 
w as a free parameter. We consider the following set of 
parameters: 

{cJb, UJc, Qs, T, Us, log[10^°As], W} . 



As stated in a previous work [18|, there exists a degen- 
eracy between the number of the extra dark radiation 
species and the dark energy equation of state, see Fig. [2] 
(upper panel). One of the main effects of a A^eff > 3.04 
comes from the change of the epoch of the radiation mat- 
ter equality, and conseque ntly , from the shift of the CMB 
acoustic peaks, see Ref. [13 for a detailed study. The 
position of the acoustic peaks is given by the so-called 
acoustic scale 9 a, which reads 




where re{zrec) and rs{zrec) are the comoving angular di- 
ameter distance to the last scattering surface and the 
sound horizon at the recombination epoch Zrec, respec- 
tively. Although ro[zrec) almost remains the same for 
different values of Aoff, rs{zrec) becomes smaller when 
A"cff is increased. Thus the positions of the acoustic 
peaks are shifted to higher multipoles (smaller angular 
scales) by increasing the value of A'off [401 ■ A dark en- 
ergy component with w > — 1 will decrease the comoving 
angular diameter distance to the last scattering surface 
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FIG. 2. The top and bottom panels depict the 68% and 95% 
CL bounds in the A'cff — w and c^;j, — w planes, respectively. 
The blue, red and green regions refer to the Baseline, Base- 
lineSPT and BaselineSPT-SNIa models, respectively. 
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FIG. 3. The top and the bottom panels depict the 68% and 
95% CL bounds in the A^eff — ns and c^g — N^h planes, respec- 
tively. The blue, red and green regions refer to the Baseline, 
BaselineSPT and BaselineSPT-SNIa models respectively. 



Teizrec), shifting the positions of the CMB aeoustic peaks 
to larger angular scales, i.e. to lower multipoles £, com- 
pensating, therefore, the effect induced by an increase of 
A^oft- The reconstructed MCMC values for w (see the 
third and fourth columns of Tab. IIIip are larger than the 
value used in the input cosmology w = —1, excluding 
the ACDM scenario with high significance. A dark ra- 
diation component which deviates from its standard be- 
havior could therefore be confused with the presence of 
a dark energy fluid with w ^ —1. 



D. TO(a)CDM model 
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We also consider models of the dark energy in which 
the equation of state of the dark energy component varies 
with time. We use a parameterization that has been 
extensively explored in the literature j4l| - |4^ : 

w(a) = Wo + Wa{l — a) , 

where wq is the equation of state par ameter at present, 
while Wa = —2dw/d\na\a=i/2 f42l. |45| . We consider the 



FIG. 4. 68% and 95% CL allowed regions in the Us — Umn 
plane from MCMC fits of Planck (blue regions) and COrE 
(red regions) CMB mock data. 



following set of parameters: 

{ujb, Wc, 6s, r, ns,log[10^°As],wo, Wa} 
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FIG. 5. Same as Fig. |4]but in the w — Us plane. 

The fifth and sixth columns of Tab. IIIII show the recon- 
structed values of wq and Wa after fitting the Plank and 
COrE mock data generated with a non standard viscos- 
ity parameter c^j^ = 0.1 but with w = — 1 to a cos- 
mology with standard dark radiation but with the time 
varying dark energy equation of state w{a) used here. 
The correlation between wq and Wa is shown in Fig. m 
The reconstructed values that we find for Planck (COrE) 
mock data are wq = —1.19 ± 0.10 and Wa — 0.77 ± 0.23 
(wo = -0.99±0.05 and Wa ^ 0.88±0.06) at 68 % CL, val- 
ues which are consistent with the current constraints on 
these two dark energy parameters, see Ref. Q- Therefore 
it is crucial to unravel the nature of the dark radiation 
component since if it turns out to be non standard, fu- 
ture cosmological data might be misinterpreted as a time 
varying dark energy fluid. 




-1.4 -1.2 -1.0 -0.8 



FIG. 6. Same as Fig. |4]but in the wo — Wa plane. 



IV. CONCLUSIONS 

Dark radiation in the Standard Model of elementary 
particles is made of three light neutrinos. However, many 



extensions of the Standard Model predict an extra dark 
radiation component parameterized in terms of the rel- 
ativistic degrees of freedom. Such is the case of sterile 
neutrinos, axions or other light degrees of freedom pro- 
duced along the thermal history of the universe. This ex- 
tra dark radiation component will be characterized not 
only by its abundance but also by its clustering prop- 
erties, as its effective sound speed and its viscosity pa- 
rameter. If dark radiation is made of sterile neutrinos it 
should have an effective sound speed clg and a viscos- 
ity parameter c^j^ such that c^g = c^jg = 1/3. However, 
other relativistic species might not behave as neutrinos, 
being c^j^ and c^g different from their canonical values. 
Current bounds on the number of relativistic species and 
on the dark radiation perturbation parameters c^g and 
c^jg have been computed using up to date cosmological 
data. We find a strong degeneracy between c^jg and the 
scalar spectral index of primordial perturbations, as well 
as between c^j^ and the dark energy equation of state w. 
The last degeneracy is alleviated when CMB data from 
the SPT experiment is added in the MCMC analyses. A 
question which naturally arises from the presence of these 
degeneracies is whether or not future CMB data will be 
able to distinguish among different dark radiation sce- 
narios. We have generated mock CMB data for the on- 
going Planck experiment and the future COrE mission 
with non standard values for the dark energy perturba- 
tion parameters. Then, we have fitted these data to a 
canonical dark radiation scenario with c^j^ = c^g = 1/3 
but with a running spectral index or with a dark energy 
component with w ^ —I, finding that non standard val- 
ues for the dark radiation perturbation parameters may 
be misinterpreted as a scale invariant power spectrum of 
primordial fluctuations or as cosmologies with a running 
spectral index or a time varying dark energy component 
with high significance. 
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